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We numerically address the issue of how the ground state topology is reflected in the finite 
temperature dynamics of the ±J Edwards- Anderson spin glass model. In this system a careful 
study of the ground state configurations allows to classify spins into two sets: solidary and non- 
solidary spins. We show that these sets quantitatively account for the dynamical heterogeneities 
found in the mean fiipping time distribution at finite low temperatures. The results highlight the 
relevance of taking into account the ground state topology in the analysis of the finite temperature 
dynamics of spin glasses. 
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Spin glass models are theparadigm of disordered sys- 
tems with slow dynamics Q The main ingredi- 
ents which define these models are quenched disorder 
and an inherent frustration in the interactions. These 
ingredients lead to a non-trivial ground state topology 
[3, and slow dynamics with spatial and dynamical het- 
erogeneities 0,[i,0,Q,SIS0|- Works analyzing the out- 
of-equilibrium properties have intuitively suggested a re- 
lation between dynamical heterogeneities and the ground 
state topology P. [ill H^. However, a quantitative un- 
derstanding of this precise relation still remains an open 
question. 

In particular, recent works have analyzed the 
dynamical heterogeneities found in three different het- 
erogeneous spin models. By studying single spin dynam- 
ics different qualitative behaviors were observed. On the 
other hand, other studies 0, H, S ^| have focused on 
spatially coarse grained quantities and analyzed hetero- 
geneities within a given coarse grained length. 

In this work we take into account a global property, 
dictated by the ground state topology, in order to analyze 
dynamical heterogeneities. We establish for the first time 
a quantitative relation between the ground state topology 
and the finite temperature dynamical properties of a spin 
glass model. We find that the dynamical heterogeneities 
are well accounted for by two sets of spins characterized 
by their role in the ground state. 

We consider in particular the two-dimensional ±J 
Edwards- Anderson (EA) spin glass model. In this model 
there exist clusters of spins which maintain their rela- 
tive orientation for all configurations of the ground state 
manifold 0, Q, . We extend this ground state in- 
formation to analyze the behavior of the system at finite 
temperatures. In order to do this we divide the system 
into two sets of spins: solidary spins, i.e. the spins that 
form these clusters, and non-solidary spins. The con- 
sequences of this division are two-fold. On one hand it 
gives us a quantitative tool to establish a relation between 



the ground state topology and the finite temperature dy- 
namical properties. On the other, it gives an intuitive 
physical frame in which to interpret the results. 

We begin our analysis considering the spin auto- 
correlation function, which clearly illustrates the different 
qualitative behaviors observed when the proposed divi- 
sion is taken into account. The non-solidary spins decor- 
relate faster than solidary spins, which suggests a rela- 
tion with the separation in fast and slow degrees of free- 
dom. In order to address this point we analyze the time- 
scale separation as observed in the mean persistence time 
and mean flipping time probability distribution functions 
0, . We show that the observed time-scale separation 
can be quantitatively accounted for by the two sets of 
spins. 

We consider the two-dimensional ±J EA model for 
spin glasses P, Q, Q, defined on a square lattice with 
periodic boundary conditions. The Hamiltonian of the 
model is 

^ = Y1 '^V^i^J (1) 

where cTj = ±1 is the spin variable and indicates 
a sum over nearest neighbors. The coupling constants 
Jij = ±J are random variables chosen from a bimodal 
distribution. The time evolution of the model is gov- 
erned by a standard Glauber Monte Carlo process with 
sequential random updates. 

In this model there exists clusters of solidary spins 
which maintain their relative orientation for all configu- 
rations of the ground state manifold 0, • This back- 
bone can be detected for each sample through the iden- 
tification of the diluted lattice J,7, 18], or its generaliza- 
tion, the rigid lattice The latter is formed by those 
bonds which are always satisfied or always frustrated in 
the ground state manifold. Notice that a backbone is 
also present in other systems such as the K-satisfiability 
model [Ulli. 
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FIG. 1: (a) A particular realization of bond disorder in an 
8x8 lattice. Single (double) lines indicate ferromagnetic (anti- 
ferromagnetic) bonds, (b) The corresponding rigid lattice 
(backbone). Full (dotted) lines indicate interactions which 
are always satisfied (frustrated) in the ground state mani- 
fold. The solidary (non-solidary) spins are indicated with 
close (open) circles. 

A particular sample of size N can be characterized by 
recognizing all its solidary spins as shown in Fig. ^ In 
order to obtain a statistical average over different real- 
izations of bond disorder, in all the results presented we 
have calculated the sets of solidary spins for 2000 differ- 
ent samples in systems with size 16 x 16. All mean values 
are obtained from averages with respect to both realiza- 
tions of bond disorder and thermal histories, as in Eq. 

121 

We begin our analysis of the out-of-equilibrium proper- 
ties by considering the two-time autocorrelation function, 

cit^^t) ^ -Y,[{MtMm (2) 

i=l 

which measures the overlap of the spin configurations at 
times tw and t . The brackets [. . .] indicate an average 
over different realizations of bond disorder, while (. . .) is 
a thermal average, i.e. an average over different initial 



conditions and realizations of the thermal noise. In each 
initial condition the spins take random values ai — ±1, 
which corresponds to a quench at t = from T = oo to 
the temperature T at which the system is analyzed. It is 
worth stressing that usually one is interested in studying 
the out-of-equilibrium properties below the critical tem- 
perature Tc- Howe ver, in the two-dimensional EA spin 
glass model Tc = [2l|, |23| ■ Nevertheless it is widely ac- 
cepted that for low enough temperatures, the dynamics 
remains out of equilibrium at short times and is very sim- 
ilar to the one observed in three-dimensional spin glasses 

[Ulli. 

For each realization of bond disorder the division in sol- 
idary and non-solidary spins can be taken into account 
rewriting the sum in Eq. ^ asC (t„, , <) — fsCs + fnsCns , 
where fg {fns — ^ ^ fs) is the fraction of solidary (non- 
solidary) spins, and Cg (Cns) is the two-time autocorre- 
lation function restricted over the solidary (non-solidary) 
spins. Note that the fraction of solidary spins is approx- 
imately 67% of the total number of spins Figure [21 
shows the behavior of C{t^,t) vs. t — when T ~ 0.5 
and tyj = 10''. For this parameters the system is in the 
aging regime |2^, and similar qualitative results are ob- 
tained for lower temperatures. The full line corresponds 
to the behavior of C when all the spins are considered. 
The behavior of Cg (Cns) is indicated with close (open) 
circles. For short times [t — < 10) the solidary spins 
are strongly correlated, i.e. they maintain their relative 
orientation in time. The non-solidary spins present a 
qualitatively different behavior, with a faster decay of the 
correlation function. Only those spins which are solidary 
in the ground state tend to remain correlated in time at 
finite temperatures. For long times {t — t^, > 10"^) each 
set of spins presents the same qualitative decay as the 
whole system. This shows that the relaxation times of 
a fraction of non-solidary spins is coupled to the ones of 
solidary spins thus decorrelating together at longer times. 
This behavior suggests a strong separation in character- 
istic times for the two sets of spins and a possible path to 
analyze dynamical heterogeneities as previously observed 
in Ref . . 

One possible path to the analysis of dynamical hetero- 
geneities is through the mean persistence time probabil- 
ity distribution function (PDF). This quantity depends 
on the time window of interest, given by i^, and t, and 
is defined as the time at which, in average, a given spin 
changes its state for the first time with respect to its 
state at tm. The mean persistence time, Tp, is obtained 
for every spin and the corresponding PDF is constructed, 
Pp{\nTp). The In Tp scale is preferred due the broadness 
of the PDF. 

In Fig. O the behavior of Pp (In Tp) is shown. The 
symbols are the same as in Fig. |21 The PDF of the 
whole system (full line) presents a sharp peak around 
In Tp 7 (Tp 10'^), with a pronounced shoulder for 
lower times. It is worth stressing that there exists a di- 
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FIG. 2: Autocorrelation function C for t^j — iC and T — 0.5. 
The full line shows the behavior of C when all spins are taken 
into account. Close (open) circles correspond to the behavior 
of the correlation Cs(C„s) when only solidary (non-solidary) 
spins are considered. 




FIG. 3: Mean persistence time PDF for the time window t — 
= 10* with = 10'' and T = 0.5. The whole distribution 
(full line) is divided in the PDF of the solidary (close circles) 
and non-solidary (open circles) spins. 



rect relation between the mean persistence time PDF and 
the autocorrelation function. For short times, the posi- 
tion of the shoulder, hi Tp ~ 2 (r^ ^ 10), corresponds 
to the first decay observed in the full C, while for long 
times, the position of the peak coincides with the second 
decay observed in the full C. The division in solidary and 
non-solidary spins allows for a physical interpretation of 
this time-scale separation. In Fig. |31we show the mean 
persistence time PDF for solidary and non-solidary spins 
separately. We observe that the shoulder found in the full 
PDF is given only by a contribution of the non-solidary 
spins. On the other hand, both solidary and non-solidary 
spins contribute to the sharp peak. The interpretation is 
straightforward: A fraction of non-solidary spins decor- 



relate first due to their low mean persistence time. At 
higher times the remaining fraction of the non-solidary 
spins and the solidary spins decorrelate together, both 
having similar mean persistence time. The same relation 
between the mean persistence time PDF and the auto- 
correlation function was observed for lower temperatures, 
giving support to our interpretation. 

We expect this particular separation in solidary and 
non-solidary spins to be reflected in other finite tempera- 
ture dynamical quantities. Recently, Ricci-Tersenghi and 
Zecchina |J| have observed a strong time-scale separation 
in the mean flipping time PDF, P/, as a signature of dy- 
namical heterogeneities. We analyze this quantity using 
the ground state information. Pf is obtained by mea- 
suring the number of flips (NfUpg) done by every spin 
within the time window extending from to t. The 
mean flipping time tj for a given tw and t is defined 
as the time window size divided by the number of flips: 

Tf = {t-t^)/NfHps 0. 

Figure 0] shows the behavior of Pf (In t/ ) . The sym- 
bols are the same as in Fig. El The PDF of the whole 
system (full line) presents two main peaks 0, which are 
a manifestation of strong dynamical heterogeneities [3| ■ 
Generally speaking these two peaks correspond to fast 
(left peak) and slow (right peak) spins. We also mea- 
sure the mean flipping time distribution for solidary and 
non-solidary spins separately. In Fig. 21we show that the 
two peaks of the full PDF can be well accounted for by 
this separation. The slow (fast) spins at finite tempera- 
ture correspond to solidary (non-solidary) spins. At high 
temperatures, the two peaks collide and the strong time 
scale separation is no longer observable. 

It is worth stressing that this separation reveals that 
a further internal structure is present, as can be seen in 
the shoulder observed in the mean flipping time PDF of 
the solidary spins in Fig. ^ For low temperatures we 
observed that the shoulder does not seem to depend on 
temperature. Instead, the peak of the slow spins moves 
to higher values in accordance with an activation process 
with a characteristic energy barrier (see inset in Fig. 0J. 
This energy barrier, 4J, corresponds to flipping a spin 
with only one frustrated bond. This should be contrasted 
with the fact that the peak of the fast spins does not move 
with temperature as shown in the inset. However, we 
must point out that a possible difference could be present 
in the tails of the distributions. For both, fast and slow 
spins, the tails seem to be power-law like, P/(t/) ~ '^^^ '^ ■ 

Summarizing, we have presented a numerical study of 
the two-dimensional ±J EA spin glass focusing on how 
the information of the topology of the ground state man- 
ifests in the finite low temperature dynamics. We have 
concentrated in the pre-asymptotic aging regime of this 
particular model as representative of glassy dynamics. 

In the ground state, spins can be divided in two sets, 
solidary and non-solidary spins. In the EA model, this 
characterization is non-trivial and deserves careful and 
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our results, and should be of relevance for understanding 
heterogeneities in EA spin glasses with continuous 
coupling distributions. 
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FIG. 4: Mean flipping time PDF for the time window t — t^ = 
10* with i„ = 10'' and T = 0.5. The whole distribution (full 
line) is divided in the PDF of the solidary (close circles) and 
non-solidary (open circles) spins. The power-law like behavior 
of the distribution's tails is highlighted. The inset shows the 
behavior of the maximum of the solidary (nonsolidary) spins 
PDF, In Tf max2 (In T/ maxi), for three different temperatures: 
T = 0.4, 0.5 and 0.6. 



time consuming simulations Once these two sets 

were identified, we analyzed the contribution of each set 
to the finite temperature dynamics. The autocorrelation 
function for each set of spin behaves differently, showing 
a faster initial decay for non-solidary spins. This natu- 
rally leads to the analysis of dynamical heterogeneities. 
First, we analyze the mean persistence time distribution, 
and show that it is intimately related to the two-step re- 
laxation of the autocorrelation function. A fraction of 
non-solidary spins, with lower mean persistence times, 
give rise to the first decay of the autocorrelation func- 
tion. The decay observed at longer times, corresponds to 
a peak in the mean persistence time distribution, and is 
shared by solidary and non-solidary spins. 

Finally, we test the relevance of the separation in 
solidary and non-solidary spins in the mean flipping 
time distribution, which presents strong dynamical 
heterogeneities. As was already pointed out in Ref. 
this distribution presents two sharp peaks. This 
time-scale separation was used to dynamically define 
groups of slow and fast spins Here we show for 

the first time that these dynamical characterization is 
well accounted for by the ground state characterization 
in solidary and non-solidary spins. Furthermore, new 
interesting and promising questions arise. For instance, 
the mean flipping time distribution for solidary spins 
presented in Fig. 0] presents a clear shoulder at low 
mean flipping times. This shoulder corresponds to an 
internal structure within the set of solidary spins. This 
suggests that a further division into sub-sets could refine 
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